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• Premise of the study: Microsatellite primers were identified and characterized in Acca sellowiana in order to expand the limited 
number of pre-existing polymorphic markers for use in population genetic studies for conservation, phylogeography, breeding, 
and domestication. 

• Methods and Results: A total of 10 polymorphic microsatellite primers were designed from clones obtained from a simple se- 
quence repeat (SSR)-enriched genomic library. The primers amplified di- and trinucleotide repeats with four to 27 alleles per 
locus. In all tested populations, the observed heterozygosity ranged from 0.269 to 1.0. 

• Conclusions: These new polymorphic SSR markers will allow future genetic studies to be denser, either for genetic structure 
characterization of natural populations or for studies involving genetic breeding and domestication process in A. sellowiana. 
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Acca sellowiana (O. Berg) Burret (Myrtaceae), commonly 
known as feijoa or pineapple guava, naturally occurs in south- 
em Brazil and Uruguay (Thorp and Bieleski, 2002). Santos 
et al. (2009) have suggested that there is potential to develop 
this fruit resource as a food product and as a medicinal plant. In 
addition to its unique taste, the species has several bioactive 
properties that make it potentially useful for pharmacological 
purposes (Weston, 2010). The species is neither well known 
nor thoroughly exploited in its center of origin; however, com- 
mercial cultivars of the species have been cultivated in other 
parts of the globe (e.g., Colombia and New Zealand). The de- 
velopment of simple sequence repeat (SSR) markers can be 
very useful for a number of studies, such as cultivar fingerprint- 
ing and genetic diversity assessment of germplasm, which are 
essential for conservation, phylogeography, breeding, and char- 
acterization of population genetic structure of A. sellowiana. 

Some specific molecular markers have already been devel- 
oped (Santos et al., 2008); however, preliminary studies using 
those markers show that only nine of the 13 markers are being 
successfully used (data not published). Therefore, the development 
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of new SSR markers will help to characterize the levels and 
distribution of genetic diversity; they will also be essential for 
phylogeographic, species domestication, and genetic mapping 
studies. Therefore, in this paper, we present 10 new microsatel- 
lite markers for A. sellowiana. 



METHODS AND RESULTS 

Total DNA was isolated from leaves of a single feijoa plant according to 
Doyle and Doyle (1990), and an SSR-enriched library was constructed based on 
Billotte et al. (1999). Briefly, total DNA was digested using Afal restriction 
enzyme (Invitrogen, Carlsbad, California, USA) at 37°C for 1 h, followed by 
enzyme inactivation at 65°C for 20 min. Digested DNA was then ligated to 
Afal adapters (Afa21: 5'-CTCTTGCTTACGCGTGGACTA-3' and Afa25: 
5'-TAGTCCACGCGTAAGCAAGAGCACA-30 in a 50-|iL reaction containing 
500 ng of digested DNA, 10 |iM of each Afal adapter, Ix T4 DNA Ligase Buf- 
fer (Promega Corporation, Madison, Wisconsin, USA), and 1 U/jiL of T4 DNA 
Ligase (Promega Corporation). The reaction was incubated at 20°C for 2 h, 
followed by enzyme inactivation at 70°C for 10 min. Fragments were PCR 
amplified using a Veriti Thermal Cycler (Life Technologies, Carlsbad, Califor- 
nia, USA) to yield adequate DNA for the next step. The reaction was performed 
in a 50-jiL reaction containing 5 jiL of the ligation reaction, 10 jiM of Afa21 
adapter, 0.2 mM of each dNTP (Fermentas, Vilnius, Lithuania), Ix Taq buffer 
(Fermentas), 1.5 mM of magnesium chloride (Fermentas), and 1 unit of Taq 
DNA Polymerase (Fermentas). A Veriti Thermal Cycler (Life Technologies) 
was used to amplify the ligated DNA with the following sequence: 95°C for 
4 min; followed by 20 cycles of 94°C for 30 s, 60°C for 1 min, 72°C for 1 min; 
and a final extension step at 72°C for 8 min. The DNA was enriched for micro- 
satellite repeat regions (CT^, GT^, and TTCJ using Streptavidin MagneSphere 
Paramagnetic Particles (Promega Corporation) using the protocol suggested by 
the manufacturer. 

Twenty microliters of the enriched DNA was PCR amplified with the same 
conditions used in the previous reaction, cloned using pGEM-T Easy Vector 
(Promega Corporation) according to the protocol suggested by the manufac- 
turer, and transformed in chemically competent Escherichia coli XL 1 -Blue 
cells (Agilent Technologies, Santa Clara, California, USA). A total of 192 posi- 
tive clones were sequenced with a MegaBACE 1000 automated sequencer (GE 
Healthcare Biosciences, Pittsburgh, Pennsylvania, USA). Of the sequenced 
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Table 1. Characteristics of 10 microsatellite primers developed in Acca sellowiana. 



Locus 




Primer sequences (5'— 3') 


Repeat motif 


Allele size (bp) 




GenBank accession no. 


Fse04 


F 
R 


<6-FAM>TTCGTTTGGCGTTTACCTTC 
CGGTGCTCGTTTGTTGTATG 


(CT)i2 


204 


59 


JX978696 


Fse06 


F 
R 


<6-FAM>CGAATAGGGCACTCCAACTC 
GCCGTCCTTAGGTTCAACAC 


(CT)9 


232 


54 


KJ14661 


Fse08 


F 
R 


<6-FAM>TCAGGTGTGGAATCTGCTTCT 
CCCTCTCATTAGGTGGTTGG 


(CT)io 


195 


53 


KF011978 


Fse09 


F 
R 


<6-FAM>ACCTTGCAAAAGGATGTCGT 
GCTGTTGGAGGAAAAAGCAG 


(TTOv 


166 


55 


JX978698 


FselO 


F 
R 


<6-FAM>TTTCGTCTCATACAAGCCTACAG 
CGAGGAATTCCCAAAACAAC 


(CT)6 


249 


55 


KJ146612 


Fsell 


F 
R 


<6-FAM>AAAAGGAAAGTACGCATCCA 
CCCAGCAACAGGTGTATGA 


(AG)i3 


236 


55 


JX978699 


Fsel2 


F 
R 


< 6 - F AM> GC G AT T G AAC CAT G AAG T C C 
CGGTGGATTCAATGGAGAAG 


(CT)8 


193 


55 


KFO 11979 


Fsel6 


F 
R 


<6-FAM>CCATTGTTTTGGAAGGAACA 
TTGCGATTTT GAGAGT GGAG 


(CT)2o 


167 


55 


JX978702 


Fsel7 


F 
R 


<6-FAM>CCATGAAATTGCTAAGCTGGA 
CGTGCTCGATGAAGAGGAA 


(CT)3o 


204 


55 


JX978703 


Fse21 


F 
R 


<6-FAM>ATCAGAAAATGCTGCCGAGT 
GCAACATTGCTCAGCAGGTA 


(CT)i2 


168 


57 


JX978700 



Note: = annealing temperature. 



clones, 36 presented unique SSR regions, and only 26 were suitable for primer 
design using Primer3 (Rozen and Skaletsky, 2000). Ten SSRs (Table 1) were 
genotyped in four populations of A. sellowiana: Ipe, Urubici, Bituruna, and 
Guarapuava (Fig. 1). The remaining designed primers were not included in this 
work due to PGR reaction problems (e.g., non-amplification, several off-target 
peaks). PGR reactions were performed in 12-pL reactions containing 20 ng of 
DNA template, 0.2 mM from each dNTP (Fermentas), 1 unit of Taq DNA Poly- 
merase (Fermentas), Ix Taq buffer (Fermentas), 2 mM of magnesium chlo- 
ride (Fermentas), and 5 pmol of each primer. A Veriti Thermal Gycler (Life 



Technologies) was used to amplify the 10 SSR loci with the following se- 
quence: 95°G for 5 min; followed by 32 cycles of 95°G for 30 s, annealing 
temperature (Table 1) for 30 s, 72°G for 30 s; and a final extension step of 7 min 
at 72°G. Genotyping was performed using the MegaBAGE 1000 automated 
sequencer, with ET-ROX 400 size standard (GE Healthcare). Sizing was car- 
ried out using Fragment Profiler software version 1.2 (GE Healthcare). MIGRO- 
GHEGKER version 2.2.4 (van Gosterhout et al., 2004) was used to assess null 
alleles and scoring errors. Deviations from Hardy- Weinberg equilibrium 
(HWE) for each of the four populations were calculated using GenAlEx 6.5 



SOUTH AMERICA 
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95 190 



380 



570 
= Kilometers 



Fig. 1. Geographic localities of the populations sampled in this paper: 1 = Ipe, Rio Grande do Sul (28°75'14''S, 51°35'39''W); 2 = Urubici, Santa 
Catarina (28°04'18"S, 50°2970"W); 3 = Bituruna, Parana (26°02'52"S, 51°70'20''W); and 4 = Guarapuava, Parana (25°26'29"S, 51°39'17''W). 
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Table 2. Results of initial primer screening in four populations of Acca sellowiana. 



Ipe, Rio Grande do Sul {n= 15) Urubici, Santa Catarina {n = 15) Bituruna, Parana {n = 15) Guarapuava, Parana {n = 15) 



Locus 


A 


He 


Ho 


PIC 


A 


He 


Ho 


PIC 


A 


He 


Ho 


PIC 


A 


He 


Ho 


PIC 


Fse04"^ 


6 


0.671 


0.667 


0.596 


9 


0.809** 


0.733 


0.761 


1 


0.744** 


0.813 


0.676 


9 


0.874** 


0.667 


0.826 


Fse06* 


10 


0.877 


0.750 


0.822 


11 


0.932 


0.800 


0.874 


5 


0.788 


0.667 


0.680 


8 


0.894** 


0.714 


0.846 


Fse08* 


5 


0.458 


0.400 


0.421 


5 


0.743 


0.733 


0.672 


4 


0.599 


0.500 


0.493 


5 


0.637** 


0.400 


0.541 


Fse09"^ 


4 


0.690** 


0.667 


0.610 


5 


0.533 


0.467 


0.479 


4 


0.651 


0.688 


0.578 


5 


0.729** 


0.467 


0.650 


FselO* 


3 


0.530 


0.333 


0.424 


2 


0.467** 


0.200 


0.332 


3 


0.654 


0.273 


0.553 


2 


0.536 


0.250 


0.359 


Fsell"^ 


7 


0.724 


0.800 


0.654 


8 


0.832** 


0.867 


0.781 


7 


0.778 


0.813 


0.720 


6 


0.699 


0.933 


0.635 


Fsel2* 


14 


0.940 


1.000 


0.902 


13 


0.906 


0.800 


0.866 


7 


0.740 


0.563 


0.686 


12 


0.887** 


0.867 


0.844 


Fsel6"^ 


12 


0.913 


1.000 


0.872 


15 


0.949 


1.000 


0.912 


12 


0.911 


1.000 


0.871 


8 


0.855 


1.000 


0.806 


Fsel7"^ 


12 


0.931 


0.933 


0.891 


15 


0.915 


0.800 


0.875 


14 


0.883 


0.875 


0.845 


11 


0.864** 


0.733 


0.821 


Fse2P^ 


7 


0.825** 


0.867 


0.769 


6 


0.848** 


0.867 


0.794 


7 


0.851** 


0.875 


0.800 


5 


0.789** 


0.533 


0.724 



Note: A = number of alleles sampled; = expected heterozygosity; = observed heterozygosity; n = number of individuals sampled; ns = nonsignificant 
values for null alleles analysis (P > 0.05); PIC = polymorphism information content. 
* Significant values for null alleles analysis (P < 0.05). 
** Indicates deviation from Hardy- Weinberg equilibrium (P < 0.05). 



(Peakall and Smouse, 2012). Number of alleles (A), expected (H^) and observed 
heterozygosities (HJ, and polymorphic information content (PIC) were analyzed 
using CERVUS 3.0.3 software (Kalinowski et al., 2007). 

All PCR products met the expected sizes based on sequence information 
from the DNA library. Considering all populations and loci, A ranged from four 
to 27 alleles per locus (average: 14.10). and H^, considering all populations 
and loci, ranged from 0.269 to 1.0 (average: 0.735) and 0.569 to 0.938 (average: 
0.805), respectively. The PIC, considering all populations and loci, ranged from 
0.495 to 0.923 (mean: 0.773) (Table 2). Out of 10 new SSR markers developed, 
four showed PIC values greater than the SSR markers developed by Santos 
et al. (2008), thus expanding the number of highly polymorphic markers for the 
species. HWE was tested in all populations. All loci were in disequilibrium in 
at least one population (except the locus Fsel6, P > 0.05), most likely due to 
heterozygote deficiency. However, these results should be treated with caution, 
given the small sample size for the populations (n= 15) (Hedrick, 2000). 

In addition, for the first time, the transferability of these 10 new SSR mark- 
ers was tested against DNA from four other Myrtaceae species — Plinia jaboti- 
caba (Veil.) Kausel, P. cauliflora (DC.) Kausel, P. trunciflora (O. Berg) Kausel, 
and Eugenia uniflora L. Six loci successfully amplified products in P. jaboti- 
caba, P. cauliflora, and P. trunciflora (Fse04, Fse09, Fsell, Fsel6, Fsel7, 
Fse21), and none of them amplified products in E. uniflora. These are the first 
SSR markers specific for A. sellowiana transferred to other Myrtaceae species. 



CONCLUSIONS 

This is the second study characterizing SSR markers for A. 
sellowiana. The new polymorphic markers double the sampling 
ability in the feijoa genome as well as the genome of other 
Myrtaceae species (due to the successful transferability of six 
markers), and will allow additional studies on breeding and the 
domestication process. Furthermore, the increase in genome 
sampling can help to better elaborate the proposition of con- 
servation measures for this important species from southern 
Brazil. 



LITERATURE CITED 

BiLLOTTE, N., P. J. L. Lagoda, A. M. RiSTERUcci, AND F. C. Baurens. 1999. 
Microsatellite-enriched libraries: Applied methodology for the devel- 
opment of SSR markers in tropical crops. Fruits 54: 277-288. 

Doyle, J. J., and J. L. Doyle. 1990. Isolation of plant DNA from fresh 
tissue. Focus (San Francisco, Calif.) 12: 13-15. 

Hedrick, P. W. 2000. Genetics of populations, 2nd ed. Jonas and Bartlett, 
Boston, Massachusetts, USA. 

Kalinowski, S. T., M. L. Taper, and T. C. Marshall. 2007. Revising 
how the computer program CERVUS accommodates genotyping er- 
ror increases success in paternity assignment. Molecular Ecology 16: 
1099-1106. 

Peakall, R., and P. E. Smouse. 2012. GenAlEx 6.5: Genetic analysis in 
Excel. Population genetic software for teaching and research — An up- 
date. Bio informatics (Oxford, England) 28: 2537-2539. 

RozEN, S., AND H. J. Skaletsky. 2000. Primer3 on the WWW for general 
users and for biologist programmers. In S. Misener and S. A. Krawetz 
[eds.]. Methods in molecular biology, vol. 132: Bioinformatics meth- 
ods and protocols, 365-386. Humana Press, Totowa, New Jersey, 
USA. 

Santos, K. L., M. O. Santos, P. R. Laborda, A. P. Souza, N. Peroni, 
AND R. O. NoDARi. 2008. Isolation and characterization of microsat- 
ellite markers in Acca sellowiana (Berg) Burret. Molecular Ecology 
Resources 8: 998-1000. 

Santos, K. L., N. Peroni, R. P. Guries, and R. O. Nodari. 2009. Traditional 
knowledge and management of feijoa (Acca sellowiana) in southern 
Brazil. Economic Botany 63: 204-214. 

Thorp, G., and R. Bieleski. 2002. Feijoas: Origins, cultivation and uses. 
HortResearch, Auckland, New Zealand. 

VAN OOSTERHOUT, C, W. F. HUTCHINSON, D. P. M. WiLLS, AND P. ShIPLEY. 

2004. MICRO-CHECKER: Software for identifying and correcting 
genotyping errors in microsatellite data. Molecular Ecology Notes 4: 
535-553. 

Weston, R. J. 2010. Bioactive products from fruit of the feijoa (Feijoa 
sellowiana, Myrtaceae): A review. Food Chemistry 121: 923-926. 



http://www.bioone.org/loi/apps 



3 of 4 



Applications in Plant Sciences 2014 2(6): 1400020 
doi:10.3732/apps. 1400020 



Klabunde et al. — Acca sellowiana microsatellites 



Appendix 1 . Information on voucher specimens for Acca sellowiana depos- 
ited in the herbarium of Universidade Federal de Santa Catarina, 
Florianopolis, Brazil (FLOR). 



Collection locality 



Geographic coordinates Voucher specimen 



Ipe, Rio Grande do Sul 
Urubici, Santa Catarina 
Bituruna, Parana 
Guarapuava, Parana 



28°75'14"S,51°35'39"W 
28°04'18"S,50°2970"W 
26°0r52"S,51°70'20"W 
25°26'29"S,51°39'17"W 



FLOR 49843 
FLOR 49841 
FLOR 49842 
FLOR 49839 
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